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Abstract: The Wilson coefficient C7 governing the radiative electromagnetic de- 
cays of B meson has been calculated to a very high accuracy in the Standard Model, 
but experimental bounds on either the magnitude or the sign of Cj are often model- 
dependent. In the present paper, we attempt at constraining both the magnitude 
and sign of Cj using a systematic approach. We consider already measured observ- 
ables like the branching ratios of B — > X s fi + fi~ and B — > X s 7, the isospin and CP 

asymmetries in B — > K*j, as well as Afb and in B — > K*£ + £~. We also discuss 

(2) 

the transverse observable A T ' which, once measured, may help to disentangle some 
of the scenarios considered. We explore the constraints on CV, Cg, C10 as well as 
their chirality-flipped counterparts. Within our framework, we find that we need to 
extend the constraints up to 1.6 a to allow for the "flipped-sign solution" of C 7 . The 
SM solution for Cj exhibits a very mild tension if New Physics is allowed in dipole 
operators only. We provide semi-numerical expressions for all these observables as 
functions of the relevant Wilson coefficients at the low scale. 
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1. Introduction 



In the last decade, one of the main avenues to search for New Physics signals in B and 
K decays has consisted in overdetermining the parameters of the Cabibbo-Kobayashi- 
Maskawa matrix (which encodes charged weak transitions in the Standard Model 
(SM)) and its representation as a unitarity triangle embedding CP- violation. The 
resulting picture has shown a very good overall agreement of all the constraints, apart 
from some discrepancies (direct CP asymmetries difference between B~ — > K~tt° 
and B° — > K~tt + , B — ^ rv versus sin 2/3, B s meson mixing from J/^cf) channel, 
and the dimuon asymmetry), which are still under experimental scrutiny but may 
be understood in terms of New Physics contributions []!], |], [|, |], |J . 

In the meanwhile, a long list of rare B decays has been determined at present 
with high theoretical and experimental accuracy. A tool of choice for these anal- 
ysis is the effective Hamiltonian describing flavour transitions, allowing an elegant 
separation between long-distance operators Oi (leading to contributions governed by 
strong and electromagnetic SM interactions) and short-distance Wilson coefficients 
Ci (summing up all the details of the fundamental theory lying beyond the SM at 
higher energies). Once expressed in this language, the analysis of rare B decays 
corresponds to constraining the allowed range of Wilson coefficients (WC), taking 
into account several observables. One must be careful that New Physics (NP) can 
not only change the value of the SM Wilson coefficients, but also introduce new 
operators with a Dirac structure that is different from the SM ones. We hope that 
overconstraining these Wilson coefficients will push them into regions incompatible 
with the Standard Model, providing hints of the structure of the underlying theory 
responsible for these New Physics effects (right-handed currents, scalar or tensor 
contributions, etc.). 

This program turns out to be quite challenging as many observables depend 
not on a single WC but a combination of many of them. Hence the constraint 
on a particular WC depends very much on the assumptions made on the type of 
New Physics present and its impact on different WCs. Many model-independent 
analyses with the aim of avoiding fine tunning assume that only the Wilson coefficient 
analysed receives a contribution from New Physics (all the other ones being set 
to their SM values). The limits of such an approach are quite obvious, and the 
conclusions that can be extracted are rather limited, specially when the framework 
is not clearly defined. As an illustration, it was proposed sometime ago to consider 
a NP contribution to the WC of the electromagnetic operator O7 approximately 
twice as large as the SM one but of the opposite sign, so that the prediction for 
B(B — > X s j) would be similar to that of the Standard Model, which was in good 
agreement with the current experimental value. This solution attracted some interest 
recently, as it could explain the Belle measurements || for the exclusive decay B — > 
K*£ + £~ suggesting that the forward-backward asymmetry did not exhibit any zero 
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at low energies. In ref. JF|, this so-called 'flipped-sign solution" was shown to be at 
odds with the prediction of B(B — > X s £ + £~). More generally, this question can be 
answered only once we fix the values of the other operators that can contribute to 
the observables: the conclusions may change if NP is allowed to contribute also to 
the semileptonic operators Og^o, or if relevant operators with a non-SM structure 
are included. Other solutions to this forward-backward asymmetry issue were also 
discussed in ref. M. 

Fortunately, the rich phenomenology of B decays together with the increasingly 
large amount of data from B factories and hadron machines open new perspectives 
to deal with larger sets of operators. In this article, we propose to focus on the two 
Wilson coefficients associated with the electromagnetic operator O7 and its chirally- 
flipped counterpart Oj> as tools to search for New Physics in a systematic approach. 
Our goal is that these coefficients play here a similar role to the p and fj parameters 
in the studies of the unitarity triangle. CV and Cy do not exhaust all the information 
that can be obtained concerning New Physics, exactly as p and fj are not sufficient 
to describe the full structure of the CKM matrix, but they provide an interesting 
summary of the situation and a good starting point to investigate NP contributions 
with other structures. 

We will focus on the allowed regions for this pair of Wilson coefficients under dif- 
ferent scenarios defined later on and corresponding to letting more and more Wilson 
coefficients receive New Physics contributions. Each scenario will be more general 
than the previous one. The basic idea is that different choices of NP scenarios may in 
principle lead to different solutions or allowed regions for each Wilson coefficient in 
agreement with all present constraints. The non-overlapping regions may be distin- 
guished thanks to additional observables, yet to be measured, providing a criterion 
to distinguish between the different NP scenarios. 

We will consider seven observables in our analysis. Six of them are believed to 
exhibit a limited sensitivity to hadronic uncertainties 1 : 

1. for inclusive decays, the branching ratios B(B — > X s ^) and B(B — >■ X s £ + £~), 

2. for B — > K*£ + £~ , the polarization fraction F^, the forward-backward asymme- 

(2) 

try Afb and the transverse asymmetry . 

3. for B — y K*^, the exclusive CP asymmetry Sk*^- This observable is not in 
the same footing of robustness as the previous observables, however its main 
theoretical uncertainties are reasonably under control. 

The list could be extended to include other future and theoretically-clean observables 
like (i = 3, 4, 5) proposed in ref. || . However for the sake of simplicity we will not 

1 Notice that even though we analysed the branching ratios for B(B — > K*^) and B(B — > 
K*£ + £~) we decided not to include them in the list, mainly due to the presence of significant 
hadronic uncertainties in form factors (see Fig. 
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include them in this paper. All of the observables above are measured with different 

(2) 

levels of accuracy except for A T \ which will be measured in the near future and can 
be used as an efficient probe to constrain the dipole operators in a different way from 
current observables. The seventh observable in our analysis, not included in this 
list, is the isospin asymmetry Aj(B — > K*^). Even though it is strongly sensitive 
to hadronic uncertainties, we include this asymmetry because of its discriminating 
power in our discussion of NP solutions. 

Our New Physics framework is defined by considering that NP enters in Oi with 
i = 7, 9, 10 (electromagnetic and semileptonic operators), together with the chirally- 
flipped operators Oy with i = 7, 9, 10. The precise definition and conventions for 
all those operators is presented in Section 2. We will split 2 this framework in three 
different scenarios corresponding to switching on NP step by step, starting from 
dipole operators and finishing with the full set of operators in the framework: 

• Scenario A. In this scenario the main New Physics contributions affect the 
electromagnetic dipole operators Gj, Or- 

• Scenario B. Here New Physics affects not only O7, Or, but also the SM-like 
semileptonic operators Og and O\o- 

• Scenario C. This is the most general case in the framework we have defined, 
where all operators 07,9,10 and Or ,9', 10' can receive NP contributions. 



This will allow us to have a better control, once confronted with data, on the im- 
pact of enlarging, step by step, the set of operators, as well as providing information 



on the effects from right-handed currents |1(J, 11], jl^, y|, [14|, [15, Tg, |17 |. Our guide- 



line in splitting the framework in scenarios will be to try to find in a systematic way 
the minimal set of operators compatible with data inside a framework (and extend 
it if necessary). Once this is done, a future step would be to find which theories can 
contribute to the selected operators. 

We will assume that NP enters only these operators, and that their Wilson co- 
efficients are real. If no solution compatible with all constraints is found at the end 
of our analysis, within our defined framework, the next step will consist in gener- 
alizing the framework to other operators (like scalars, tensors, the chromomagnetic 
operator 3 or further chirally-flipped operators). The generalization is systematic and 
straightforward and will be presented elsewhere, but some details will be given here. 
We classify our observables in three categories: 



2 This splitting is not unique and different choices are possible. The only condition is to start 
from a restrictive NP scenario, where only dipole operators are affected by NP, and end up with 
the most general scenario. 

3 This generalization maybe particularly interesting because it would affect most of the observ- 
ables described here, and only weak bounds on this operator are available till now. 
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1. Class-I observables mainly sensitive to 7 and O?/, but not to C^=9,io,9',io'- 

2. Class-II observables exclusively sensitive to 7 and 7 >, to semileptonic opera- 
tors (Oq and Oiq) and their chiral counterparts O w >. Only these operators 
intervene, even within more general frameworks than the one considered here. 

3. Class-Ill observables that are also sensitive to all the previous operators Oi 
with % = 7 ... 10', and in addition have the potential of exhibiting a sensitivity 
to NP contributions from other operators like scalars, tensors, chromomagnetic 
operator... 4 including all the previous operators Oi with i — 7, 7', 9, 9', 10, 10' 
but also scalar, tensor, chromomagnetic, etc., operators. 

Notice that, strictly speaking, within our defined framework, Class-II and Class-Ill 
observables coincide. However, having in mind a systematic forthcoming generaliza- 
tion of this work we need to split them as a function of their potential NP sensitivity 
beyond our present framework. We will also discuss how this classification would 
change if we extend the framework to include additional Dirac structures. It would 
basically require to re-classify some observables (mostly in Class I) and introduce 
more Class II subdivisions (even if not required here, one could also add intermedi- 
ate stages between Class II and Class III at will). 

In this paper we will illustrate the method on the practical example of determin- 
ing the sign of C 7 , already discussed in ref. JFj, using a subset of our observables 5 . 
We will focus not only on the restrictive "flipped-sign" solution, but allow also for 
deviations in the modulus of C 7 . As it is well known, the sign of C 7 has an impor- 
tant impact on observables like the forward-backward asymmetry {Apb) in the rare 
exclusive semileptonic decay B — > K*fi + fi~, that is at present slightly at odds with 
the SM prediction. 

In section 2, we present in detail the operators entering our framework and the 
observables of interest, with their current experimental accuracy as well as numerical 
expressions for the implementation of their theoretical determination. In section 3, 
we discuss the three different scenarios and combine the present constraints for each 
of those scenarios to look for different solutions or allowed regions in the WC planes. 
In section 4, we summarize the elements learned concerning the sign of C 7 and the 
values of the WCs. Most technical details concerning the inputs and the computation 
of the observables are collected in the appendices. 

4 Thcre is an important distinction in our way of treating Class-I observables with respect to the 
other classes: the definition of Class-I observables involves only their sensitivity to dipole operators 
and their lack of contributions from semileptonic operators. Other potential sensitivities beyond 
the defined framework are not relevant at this stage. This is essential to be able to define primary 
regions in a systematic way for each framework. On the contrary, we prefer to split Class II 
from Class III, to identify more easily the observables that will change if new sources beyond the 
framework are included. 

5 An interesting analysis was also presented in ref. p!q |, considering another subset of our observ- 
ables but adding NP to one Wilson coefficient at a time. See Sect. 4 for further details. 
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2. Operators, method and observables 



2.1 b — > s effective Hamiltonian 



We consider the effective Hamiltonian for radiative b — > s transitions [|19|, 



V2 

(s) 

with the CKM matrix combinations X q = V qb V* and 

6 10 



H§ = G x O\ + C 2 O c 2 + CiOi + + c i'°i')> 

i=3 i=7 

?4ff = C x (0{ - Oi) + C 2 (O c 2 - O u 2 ). (2.2) 



= C^t) (fib) and = O^o (fit,) are the Wilson coefficients and the local effec- 
tive operators respectively. The contribution of is usually dropped for being 
doubly Cabibbo-suppressed with respect to that of H^g, but we will keep it for the 
observables of interest. In eq. ( |2.ip we use the same operator basis as ref. pi| . We 
focus our attention on the operators 

7 = -^—m b (sa, u P R b)F^, O r = -^—m b (sa, u P L b)F^, 

2 2 

°* = 4^(^Pifi) (h^) , o v = -?—(s lfl p R b) (hn) , 

2 2 

°™ = ^(si.PMhV), Ow = ^— 2 (si»P R b)(h»ld), (2.3) 

where Pl,r — (1 T 7s)/2 and m b = m b (n b ) denotes the running b quark mass in 
the MS scheme. The primed operators, with flipped chirality with respect to the 
unprimed ones, are either highly suppressed or vanish in the SM. Hence, 

c sm = ^C 7 SM , Cl,% = (2.4) 
In the following, we will assume that only the Wilson coefficients of the operators in 



eq. (|2.3| ) are potentially affected by NP according to our framework. 



The determination of the Wilson coefficients in the Standard Model follows the 



discussion in refs. |19], |2(| to perform the matching at the high scale /Jo (potentially 
affected by short- distance NP) and the running of the Wilson coefficients from the 
high-scale down to fj, b , leading to SM Wilson coefficients at NNLO accuracy. The 
error budget of the observables includes a variation of fi b from twice to half its 
central value (we take fi b = 4.8 GeV). We have also checked that the variation 
of the high scale /xq yields only a tiny uncertainty on the observables. We follow 
refs. E5| and include QED corrections through five additional operators 
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cm 


C 2 M 


C 3 (M6) 




C 5 (Mb) 








Cg(Mb) 


CloOb) 


-0.2632 


1.0111 


-0.0055 


-0.0806 


0.0004 


0.0009 


-0.2923 


-0.1663 


4.0749 


-4.3085 



Table 1: NNLO Wilson coefficients in the Standard Model at the scale /x&=4.8 GeV, 
obtained from the inputs in table ||. For the computation of the observables, we 
considered a variation of /x& from half to twice its value. 



(^3,4,5,6(3 an d Ob) mixing with the ones displayed in eq. ( |2.1| ). The values of the 
Wilson coefficients at the low-scale /x& = 4.8 GeV are given in table |], where the 



definitions |26] 



r ,..„ r 1 r 4 20 80 

^8 + C 3 -^C 4 + 20C 5 -^C 6 
o 3 



Yeff 



have been used, since CV and C& always appear in these particular combinations with 
other Cj in matrix elements. 

In tables [L] and |2| we present the most important inputs used in our observables 
including the values of the Wilson coefficients in the SM. 



Mh = 4.8 GeV 


Mo = 2M W fig 


m B = 5.27950 GeV 


27 




m K * = 0.89594 GeV M 


m Bs = 5.3663 GeV 


27 




to m = 0.105658367 GeV |2J 


sin 2 W = 0.2313 

M w = 80.399 ± 0.023 GeV 


27 
27 




M z = 91.1876 GeV §7j 


a em (M z ) = 1/128.940 


19 




a s (M z ) = 0.1184 ± 0.0007 §3] 


m poio = 173 3 ± x.l GeV 
mf s {m c ) = 1.27 ± 0.09 GeV 


28 
27 




TOj S = 4.68 ± 0.03 GeV foj 
mY^{2 GeV) = 0.101 ± 0.029 GeV |tJ 


Xckm = 0.22543 ±0.0008 
p = 0.144 ±0.025 


29 
29 




4ckm = 0.805 ± 0.020 |29j 
fj = 0.342 ±0.016 || 


B(B ->• X c eD) = 0.1061 ± 0.00017 
A 2 = 0.12 GeV 2 


19 
19 




C7 = 0.58 ±0.016 Q 


A h = 0.5 GeV 

f K * || = 0.220 ± 0.005 GeV 

a(0) = 0.3ll°; 2 ° 

01,11,1.(2 GeV) = 0.03 ±0.03 

Ab(m7() = 0-51 ± 0.12 GeV 


31 
30 
64 
30 
30 




/b = 0.200 ± 0.025 GeV |0| 
/k-, ±(2 GeV) = 0.163 ±0.008 GeV Q 
^11 (0) = 0.10 ±0.03 pi 
a 2,||,j_(2 GeV) = 0.08 ± 0.06 @ 


f B , = 0.2358 ± 0.0089 GeV 


29 




t b , = 1.472 ± 0.026 ps |2j 



Table 2: Input parameters, based on refs. []19|, p7|, [E3], p9|, p0|, [pi 



2.2 Method 

We start by describing in full detail how the method applies to our previously defined 
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framework (New Physics allowed only in electromagnetic dipole, semileptonic oper- 
ators, with SM and flipped-chirality structures). We will proceed in the following 

way: 

1. We start by classifying observables in three classes, as already mentioned: 
Class-I observables sensitive only to Oj and Or contributions, Class-II ob- 
servables exclusively sensitive to the full set of operators that we consider may 
be affected by New Physics (O7, Og, O\o as well as their flipped chirality coun- 
terparts) and Class-Ill observables, not only sensitive to all these operators, 
but also to further new operators (scalars, tensors, etc). 

2. We define a reference frame of allowed regions using observables sensitive to NP 
only through a pair of Wilson coefficients, in our case C7, Cji. These reference 
regions, that we will call primary regions, are determined by Class-I observ- 
ables, and are the maximally allowed regions inside our defined framework. 
They can only shrink when new observables are added. In principle, one could 
define a different reference frame, where NP enters only in C7, but this can be 
inferred directly from the projection of our reference region along the Cj axis. 
We are in Scenario A. 

3. We add a larger set of observables with sensitivity to larger sets of operators 
(Class II and Class III), but still inside Scenario A. Those new observables 
when restricted to the (C7, Cr) plane may cut further the primary regions, 
defining a smaller allowed region inside the primary ones. 

4. In order to expand again these allowed regions (with the maximal area always 
defined as the primary regions), we will now move to Scenario B and C, allowing 
new contributions for the extra coefficients, in our case, C 9 , C w and the flipped- 
chirality ones. 

The same procedure should be repeated if other structures are included defining 
an extended framework. A discussion can be found in Sec. |3.6| . 

We will list the observables of interest for our analysis, providing in each case 
a semi-numerical expression for the observables with their central values and their 
uncertainties in the Standard Model, as well as their dependence on the deviation 
5Ci = Ci — Cf M at the low scale /i^. This treatment assumes that the analysis of 
uncertainties performed in the SM is not significantly affected by the presence of NP. 

In many places along this paper we will refer to the correlation between pairs 
of WCs, sometimes denoted as (Ci,Cj) or (5Ci,5Cj). The relation between both is 
linear Cj(/ib) = Cf M ((J>b) + $Ci- In all cases we will plot only the correlation between 
(SCi, SCj). 
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0(0,0) = 3.15 5 a = 0.23 


a (0i7) = -14.81 


<2( 7 j) = 16.68 


O(o,7') = —0.23 



Table 3: Coefficients describing the dependence of B{B — >■ X s ^) on C 7 and C 7 i. 



2.3 Class-I observables 



Class-I observables receive contributions from O7, Or but not from the semileptonic 
operators 09,10 or CV,io'- Three observables considered here fall into this category: 
the branching ratio of the inclusive radiative decay B — > X s j, as well as the isospin 
asymmetry (Aj) and the CP-asymmetry (Sk*~/) of the exclusive decay B — > K*^. 

• B[B — > Xg'y) is one of the cleanest observables in B physics from the theoret- 
ical point of view. Apart from contributions to the chromomagnetic operator, it is 
only sensitive to electromagnetic dipole operators, without pollution from other New 
Physics contributions. The currently available experimental world average is 



B(B -> X sl ) 



B 7 >1.6GcV 



(3.55 ±0.24 ±0.09) x 10" 



(2.5) 



The following formula updates the expression in ref. [20], using ref. p4] |, based on 
the NNLO SM results of P], 135], 1361 (more details can be found in App. H). 



B(B -> X sl ) 



£ 7 >1.6 GoV 



a (0 , ) ± $a ± 0(7,7) [(SC 7 ) 2 ± (6C r ) 2 ] ± 



±a(o,7) $C 7 ± a(o,7') 



10 



-4 



where the scale of the NP contributions to the Wilson coefficients SCi(fib) (with 
i = 7, 7') is taken at /z^ = 4.8 GeV. The coefficients Oj are collected in table [3|, from 
which one can extract the SM prediction, in good agreement with the experimental 
measurements: 

B(B -> X s7 )^ >16GeV = (3.15 ± 0.23) • 10" 4 (2.6) 

• Aj(B — > K*j): The measurement of the isospin asymmetry (Aj) in B — > K*~f 
was reported by BaBar and Belle, with a slightly larger neutral decay rate and hence 
a positive Ai. 



Aj = 



T(B° -> K*°7) - r(S" -> IP-7) 



(I -R+'°t + /t q -l)/2 (2.7) 



r( J B° -> k*°7) ± r(B- ->■ if*-7) 

where the two isospin-breaking ratios are / = i3(i?° — > K* 0/ -f)/B(B~ — > K*~ / y) and 
R+/o = r(T(4s) B+B-)/r(T(4s) 5°5°) . 

The recent update of BaBar collaboration |37]] with a five times larger sam- 
ple than their previous result has moved substantially Aj in the positive direction 
(A 1 = 0.066 ± 0.021 ± 0.022), being now consistent with zero at more than 2 a, while 
previously the consistency was below 1 a. The older result from the Belle collabo- 
ration PSJ: Aj = 0.012 ± 0.044 ± 0.026 requires an update to determine whether it 
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c = 4.11% 


8c = 2.52% 




d^ — —2 51757 


e (0,0) = 1 
e(o,i) = -0.0919061 
e ( o,2) = 7.49847 


e (i,o) — —5.0165 
e (2 , ) = 6.30856 



Table 4: Coefficients describing the dependence of Aj (B — > K*j) on C 7 and C 7 r. 



follows the same trend as BaBar. The average of these two measurements according 



to the Heavy Flavor Averaging Group is ||32|| : 



Af p (B -»■ iT 7 ) = 0.052 ± 0.026 . (2.8) 

In the Standard Model, Aj vanishes in naive factorisation, and it gets contribu- 
tion only from non-factorizable graphs where a photon is radiated from the spectator 
quark. This quantity was first calculated in the SM within the QCD Factorisation 
(QCDF) framework in ref. [31| and confirmed in ref. [j39fl , with a result 9.3I32 % 



Later on, it was reevaluated adding some (Cabibbo-suppressed) annihilation contri- 
butions but changing the factorisation scale from around 2 GeV to near 4.8 GeV, 
due to the fact that, below this scale, the four-quark operators factorise, so that the 
gluon exchange responsible for the running of these operators does not probe small 
scales and thus does not induce running below fit P0|. 



This observable is dominated by l/m& corrections inducing important hadronic 
uncertainties, but we include it because of its particular sensitivity to CV and CV 
which will prove very important in our discussions. The corresponding numerical 
expression is: 



MB -> iT 7 ) = c x - ^ k 7£~ 7 > v ± Sc . (2.9) 

where the non-zero coefficients are collected in table |], out of which one extracts the 
SM prediction: 

A^B ->■ K*j) SM = 0.041 ± 0.025 (2.10) 

once again in good agreement with the experimental value. 

• Sjc*^' The radiative decay b — > 57 constitutes a major probe of both the flavour 
structure of the SM and NP. In the SM, the left-handed structure of the weak inter- 
actions makes the emitted photon mainly left-handed in b decays and right-handed in 
b decays, as can be seen from the structure of the (dominant) electromagnetic dipole 
operator sumcr^bRtL)- The needed helicity flip of one of the external quarks results 
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into a factor rrib for &r — > Sl^l and a factor m s for bi — > sr^r. Therefore, at LO in 
the SM, the emission of right-handed photons is suppressed by a factor m s /mb. This 
suppression can be overridden in a large number of NP scenarios where the helicity 
flip occurs on a internal line, which may cause appearance of a factor much larger 
than m s /mb- 

The photon helicity is difficult to probe directly, but can be accessed indirectly 
using the time-dependent CP asymmetry in B° — >• K 



*o 



7: 



A 



CP 



T(B°(t) -» K*°>y) - T(B°(t) -> ir*° 7 ) 
r(Z?°(t) ->• 7?*°7) + r(5°(t) ir*° 7 ) 



5i<-* 7 sin(AmBt) — Ck*j cos(Amet), 

(2.11) 

where and K*° are observed through their decay into the CP eigenstate KgTT 
and B° mixing is assumed to be SM-like 6 . The helicity suppresion of right-handed 
photons make Aqp dominated by 5-meson mixing in the SM, irrespective of hadronic 
uncertainties. Since NP can relieve this suppression, eq. ( [2.1 1|) is a good candidate 



for null-tests of the SM [^T| [|2|, |43], . In the present article, we will focus on Sk*^ 
in eq. fl2.HI) , as it involves the interference of photons with different polarisation and 
provide interesting constraints on Cy (see Appendix |B.3| for further details). 

The experimental results available from the B factories for Sk*-j are the following: 



riexp 



-0.32^;1 (stat.) ± 0.05 (syst.) Belle |5| (535 • 10 6 BB pairs), 
-0.03 ± 0.29 (stat.) ± 0.03 (syst.) BaBar || (467 • 10 6 BB pairs), 



with the HFAG average |2| 



= -0.16 ±0.22. 
A numerical expression for this observable with our inputs is: 



5" 



Sk*~/ — f A + 



(2.12) 



(2.13) 



where / corresponds to the SM central value and 8j, Sj the corresponding error bars. 
The non-vanishing g and h coefficients can be found in table [5[ One can see that the 
SM prediction is: 

Sf A / 7 = -0.03 ± 0.01 (2.14) 



2.4 Class II 

In this set, we find some of the observables constructed out of the coefficients of the 
angular distribution of B — > K*(—t Kn)£ + £~ for which the hadronic uncertainties 

6 This assumption is compatible with the latest measurements of the CP-violating parameter 
\p/q\ — 1.0024 ± 0.0023 [[32| derived from the data gathered at B factories only. 
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-0.0297336 




0089893 


3(0,1) 


— 4-1 ^9 77zL 




— QQQQ 


3(i,o) 


= —3.17764 


^(0,1) 


= —4.51218 


3(i,i) 


= -415.441 




= -214.866 


3(0,2) 


= +8.63917 


h(0,2) 


= +290.553 


3(2,0) 


= +8.63917 


h(2,0) 


= +290.553 



Table 5: Coefficients describing the dependence of Sk^ on C7 and Cy. 



due to form factors cancel largely, and which are only dependent on some of the spin 
amplitudes involved in this decay. The observables called (with 2 = 2,3, 4, 5) fall 
inside this category 

• . This is the only observable which has not been measured yet and is 



included in our analysis though. Its unique sensitivity to CV,9',io' (shown in [0, 



|4"7j ]) and the very limited hadronic uncertainties attached to it makes it into a very 
appealing observable to distinguish between different NP scenarios. 
Its definition in terms of spin amplitudes is [2T|: 



^ 2 >=Sw (2 ' 15) 

where A± and A\\ are the corresponding spin amplitudes of the K* and q 2 (or s 
in the following) is the lepton-pair invariant mass squared. This asymmetry avoids 
one of the main sources of uncertainty for observables based on the B — > K*£ + £~ 
decay, namely the soft form factors £j_ and £n [0- is constructed to cancel its 
dependence on £_l(<? 2 ) exactly at LO and displays only a very mild sensitivity on it 
at NLO in QCDF. Its extraction from the uniangular distributions is described in 

Appendix |B.5.3f . 

(2) 

has been computed in QCDF at NLO using our inputs in table 2, the 



soft form factors described in Appendix |B.5.2j and an estimate of A/mb suppressed 



corrections of order 10%. A detailed discussion of its sensitivity to some of the 
operators in our framework can be found in ref. 0. 

After computing this asymmetry with our inputs, we have fitted the results to a 



7 The other asymmetries (with i = 2,3, 4, 5) require the determination of the full distribution 
(A$p is particularly sensitive to O10', whereas A^' 4 ^ probe the longitudinal spin amplitude). 
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Table 6: Coefficients of the polynomial functions -F(o,o) an d ^(0,0) entering SM predic- 
tion of and those of the polynomials Pi and P 2 corresponding to the associated 
upper and lower error bands respectively The second row in this table and the 
following ones indicates the dimension of the coefficients in each column. 
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Table 7: Coefficients of the polynomial functions entering A ( T ' . 



simple parametrisation of the following form 

4V) = < C V)!£;S (2.16) 
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(2) 

Table 8: Coefficients of the polynomial functions Guj^ entering . 



with the central value 



A 



(2),CV ( 2 



Ei=0,7, 



7', 9,9', 10,10' Z^jr=i,..10 



Ei=0,7,7',9,9', 10,10' J2j=i,.. 10' G(ij)(q 2 )SCiSCj 



(2.17) 



where we have introduced the definition 5Cq = 1 to write down the constant and 
linear terms in the same way as the quadratic ones. The errors on the asymmetry 
are given with respect to the SM central value F(o t o)/G(p,oy- 



X I 2\ _ D / 2\ F (0,0)(q ) 

G (o,o)(r) 



F {m (q 2 ) 

C(o,o)(g 2 ; 



P2{q 2 )- 



(2.18) 
(2.19) 



All the above functions of q 2 = s have been fitted to polynomials in this variable. 
The coefficients corresponding to the functions -^(0,0), ^(0,0), -Pi and P 2 are given in 
table |] and that of F^j) and G^j) in tables [7] and || respectively. All these coefficients 
are dimensionful (but can be easily turned into dimensionless quantities once F, G 
are expressed as a function of s = s/m 2 B ) with the dimension indicated in the second 
row of table 6 s . 

All entries of the matrices F and G should be taken to be zero, except for 
those provided in tables || [7| and |^ and those related to them through the following 
equations 



F(T,T) 



(7,7), 



F(T,9') 



'(7,9), 



F (7',9) 



F(7,9') and F( 9 / )9 /j 



-^(10,10') : 
P(10',10') 



'(9,9'), 
P(10,10) 



(9,9)- 



(2.20) 



3 As an example of how to read those tables, we provide here the function i 7 (o j o)(s): 



^(0,0) (s) = + 12904.2 - 17256.7 GeV~ 2 xs + 10543.8GeV~ 4 x s 2 - 3519.19GeV~ 6 x s 3 
+ 667.247 GeV~ 8 x s 4 - 67.3536 GeV~ 10 x s 5 + 2.78209 GeV" 12 x s 6 . 
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Most of these symmetries, and the following ones between different Fuj) ((2(ij)) 
elements, are easily understood once the large recoil limit of the spin amplitudes is 
inserted into the definition of the observable |HJ] (see Appendix |B.5.4| for details). 
Similarly for the Gr(jj) functions we have 

C(7',9) = £(7,9')) G(10,10') = <2(9,9'), 

(2(7', 7') = (2(7,7), Gfjftf) = G(7,9) an d C(9',9') — (2(10', 10') = G(10,10) = (2(9,9) (2-21) 

together with the relations between the G^j) and F^ functions: 

(2(0,7') = F(pj), (2(7,9') = F(j,9), (2(0,9') = ^(0,9) > (2(0,10') = ^(0,10) > 

(2(0,7) = F(oj'), G(7,9) = ^(7,9')j (2(0,9) = ^(0,9')) (2(0,10) = ^(0,10')- (2.22) 

In conclusion, the total number of non-zero entries of the matrices i^ij) and G^j) 
entering eq. ( |2.17| ) is 20 for each matrix 9 . 



As stated earlier, there is no current measurement of this asymmetry, but we 
will present in Sec. || the predicted value for each of the allowed regions in our 
different scenarios. 

2.5 Class III 

Here we consider observables affected by 7 , Of, Og iW , Cg^io', and in principle other 
kinds of NP operators such as scalars or tensors. The most important observable in 
this category is B{B — > X s £ + £~~) due to its limited sensitivity to non-perturbative 
physics. In the same category fall also other observables defined through the angular 
distribution of B — > K*(— > K7t)£ + £~, in particular the forward-backward asymmetry 
Ifb and the longitudinal polarisation fraction F^. 

• B(B — > X s {i + fi~) will be used only in the low-g 2 region (from 1 GeV 2 to 6 
GeV 2 ) as the theoretical prediction in the high-g 2 (above 14.4 GeV 2 ) region suffers 
from further theoretical uncertainties E9L pOl pi], SI- In the low-g 2 region, the 



branching ratio is measured to be [52 : 



(1.49 ± 0.50±°oil) x 1(T 6 (Belle) , 
B(B X s /i + fM~)i ow . q 2 = <j (1.8 ± 0.7 ± 0.5) x 10~ 6 (BaBar) , (2.23) 

1.60 ± 0.50) x 10~ 6 (Average) . 



The SM prediction for B(B X s /j, + fj,~) is (1.59 ± 0.11) x 10~ 6 |2|. With our 
inputs and including m s -suppressed terms (see App. |B.4| for more details), we obtain 



9 For instance, the 20 non-zero elements for the matrix Fuj\ correspond to the values of 

= { (0, 0), (0, 7), (0, 7'), (0, 9), (0, 9'), (0, 10), (0, 10'), (7, 7), (7, 7'), (7, 9), (7, 9'), (9, 9), (9, 9'), 
(7', 7'), (7', 9'), (7', 9), (9', 9'), (10, 10), (10, 10'), (10', 10') } 
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Table 9: Coefficients describing the dependence of B{B — > X s fi + fi ) on CV^io and 

CV',9',10'- 



the corresponding expression for the integrated branching ratio at the scale ^ 
GeV in the low-g 2 region (from 1 to 6 GeV 2 ): 



4.8 



B(B ->• X s h + ij,-) = 10 



-7 



E 

i,i=0,7,7',9,9', 10,10' 



b^j^SCiSCj ± 5 b 



The values of the non- vanishing coefficients b are listed in table 
• Ap^(q 2 ). The forward-backward asymmetry in — > K*°i 



dT/dq 2 



d(cos9[) 



d 2 T 



dq 2 dcos6i 



d 2 T 
_! dq z dcost>i 



(2.24) 

is defined by: 
(2.25) 



with Q\ the angle between the positively charged lepton in dimuon rest frame and the 
direction of the dilepton in the Bd rest frame. This asymmetry can also be written 
in terms of spin amplitudes j|7| inside our framework as 10 



-4 FB (g 2 



r" 1 dT/dq 2 



[Re{A\\ L A\ L ) - Re(A llR A* ±R 



(2.26) 



(See Appendix |B.5| for definitions). The overall minus sign with respect to eq. (4.4) in 
ref. ^7] stems from the definition of Ap-Q(q 2 ) in eq. ( |2.25[ ) chosen to match the plots in 
refs. || G3J. The expression of dT/dq 2 in terms of K* spin amplitudes (including the 

muon mass terms) can be found in eq. (|B.42|) . The QCDF framework at NLO is well 

(2) 

suited to compute A FB , just as we did previously with AL , including an estimate of 
A /mj corrections. However, unlike A\, , Afb can receive not only contributions from 

the operators Oi, 0[ with i — 7, 9, 10 but also from scalar and tensor operators || El]. 

(2) 

Another important difference between A^ and A-p-Q is that Afb is not protected at 



LO from soft form factor uncertainties contrary to Ajp . Besides, A^' exhibits the 



I (2) 



10? 



^Notice that while eq. ( 2.25| ) is valid in general, eq. (2.26) is valid only within our framework, 
which means that one should add extra amplitudes in eq. (2.26) when scalar operators are included. 
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Table 10: Coefficients of the polynomial functions #(o,o) an d -f(o,o) entering SM predic- 
tion of A FB and those of the polynomials P 3 and P 4 corresponding to the associated 
upper and lower error bands respectively 



same remarkable features as Aps like, for instance, the presence or absence of a zero 
(in the presence of right-handed currents) || ^BJ. A FB has been under scrutiny 
lately, as a consequence of the Belle measurement suggesting that, contrary to SM 
prediction, it might not display a zero in the low-g 2 region, triggering many proposals 
to explain this behaviour [|[ Q. 

We define the integrated forward-backward asymmetry in the low-g 2 region to 
agree with the experimental determination: 

A " B ~ ,6&,V jr ' (2 - 27) 

JlGeV 2 <2g 2 



while the average of the measured values by Belle || and CDF collaborations |53j is 

A e *» = 0.33^; 2 i (2.28) 



We can provide a semi-numerical expression for this observable in a similar way 
to A T . Starting from the unintegrated asymmetry 

A FB {q 2 ) = A°%{<?t\%], (2.29) 
where the central value (CV) is 

AGV ( 2\ _ ^i=0,7,7'.9,9'.10,10' ^2 j=i,..., 10' H(ij)(q 2 )SCjSCj 
FB l<? ) — "y^ 7 / 2 ur< 77^ (/.oUJ 

2^i=0,7,7',9,9',10,10' l^j=i,...,10> I (hj) W )°^i ^j 

(using again 5C*o = 1) and the uncertainties are given with respect to the SM central 
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Table 11: Coefficients of the polynomial functions Huj\ entering A-p-g. 



value curve (if( 0i o)//(Q,o)] 



- Pa(g 2 ) - (2.31) 
W)-f^-^ 2 )- (2-32) 

After integrating over the low-g 2 experimental kinematic range (1 < q 2 < 6GeV 2 ), 
following eq. (|2.27|) we obtain 



Afb = A c f v b % (2-33) 
where the central value can be split into SM and NP contributions: 

A% = A™ + A%, (2.34) 

with 

asm _ JiGeV* H m (q 2 )k(q 2 )dq 2 

FB ~ r 6GeV 2 ' (^-00) 

JlGeV 2 J (0,0)(? 2 )^(9 2 )^ 2 

*np _ iiGeV 2 zCi=o,7,7 , ,9,9 / , lo.io' zCj=i,..io' H^^(q 2 )k(q 2 )5Ci5Cjdq 2 { 

■"■FB — 6GeV 2 A FB) 1^.<J0J 

JiGeV 2 zCi=o,7,7',9,9', 10,10' zCj=j,.. io' hi,i){<l 2 )k{q 2 )bCi5Cjdq 2 
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Table 12: Coefficients of the polynomials functions entering ^4fb- 



and the uncertainties are defined, according to eq. fl2.32| ), as 

° u ~ ^75^ ......... ' V- 6( ) 



JiGev 2 hmi^M^W 2 
Iigow 2 7 (o,o)(g 2 )% 2 )dg 2 



(2.38) 



The coefficients of the polynomials #(o,o), -^(0,0), Pa an d -P4 can be found in table |10 
and those of Hn^ and are in tables [XT] and [12| respectively. 11 



11 Notice that a normalization factor fc(g 2 ) has been introduced in Eqs.( 2.35 )-( 2.38 ). This factor 
cancels exactly in the differential ApB(q 2 ) or F^(q 2 ) expressions, but should be considered for 
the integrated A FB and F L observables. Its value is: k(s) = 64.7723 - 76.8789 s + 45.3645 s 2 - 
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All components of the matrices H and / are taken to be zero (as it was done 

(2) 

for Aj. ) except for those provided in these tables and those related to them via the 
equations 

#(7',10') — -#(7,10) > #(9',10') = -#(9,10) an d P(7',10) = -#(7,10')- ( 2 -39) 

and 

I(V,W) — ^(7,9), -f(7',9) = ^(7,9')> 

^(7',7') = ^(7,7)j ^(9',9') = ^(9,9) an d PlO',10') = Pl0,10), (2.40) 



which leaves finally 12 Huj\ and 20 Iuj) non-zero functions entering eq. ( 2^30 ) 
eq. (pi). 



Using eqs. (|2.35|) , (|2.37|) , (|2.38|) and table [H] we get the following prediction for 



the integrated forward-backward asymmetry (Afb) i n the SM: 

A s F M s = -0.0316lg.-gS. (2-41) 

• Pl: The longitudinal polarization fraction of the K* in the exclusive B — > 
K*£ + £~ decay is defined in terms of the spin amplitudes as 

(2.42) 



IA) 


2 


or 




dq- 




54 





can also be computed in QCDF and, as before, an estimate of A/mb corrections has 
been added to the other sources of uncertainty of this observable. 

The integrated version of this observable in the low-g 2 region can be defined as 
in eq. (ggTfl 



r 6GcV 2 dr 



(2.43) 



rOLre V dT_ 
JlGeV 2 dq 2 

and the average of the data measured by Belle M and CDF collaborations |k| from 
this observable yields 

pf p = 0.60i°;^ . (2.44) 
The analysis of A FB and A FB performed in eqs. (|2.29|) - (|2.38|) can be repeated, 



step by step, for Pl and Pl with the substitutions P(o,o) — * J(o,o), #(i,j) — ► J{i,j)i 
P3 — > P5, P4 — > Pq and, obviously, A-pB — > Pl, ^4fb — > Pl and exactly the same 
normalization factor k(q 2 ). Table [13] contains the coefficients of J(o,o), I(o,o) (for 



11.9407 s 3 - 0.970171 s 4 + 1.79656 s 5 - 0.631163 s 6 + 0.120215 a 7 - 0.0135664 s s + 0.000856554 s 9 - 
0.0000234364 s 10 with s = q 2 . This normalization function is obtained by requesting that the SM 
prediction for the differential decay rate is exactly reproduced, shown in Fig.|l5|, i.e., 10 21 x dT /dq 2 = 
Ao,o)(9 2 )xfc(q 2 )xl0- 6 . 
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1 
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Lie v 


Li6 V 


Li( V 


\jB V 


vje V 


"•'(0,0) 


^ O \J . 1 




1^ ± O i O -LO 








+ 1 1 7 1 02 


7 (0,0) 


+773134 


-72762.1 


+280788 


-88514.3 


+24423.2 


-3375.38 


+ 188.857 


p 5 


-.0792139 


+.952685 


-.395205 


+.0821238 


-.00911051 


+4.67994- 10- 4 


-6.09404- 10" 6 


p 6 


-.133068 


+.720264 


-.154064 


-.0186277 


+.0121348 


-1.77815- lO" 3 


+8.87194- 10" 5 



Table 13: Coefficients of the polynomial functions «/(o,o) an d /(o,o) entering SM pre- 
diction of Fl and those of the polynomials P5 and P$ corresponding to the associated 
upper and lower error bands respectively 



(U) 


1 


s 


s 2 


s 3 


s 4 


s 5 


s 6 


dim 


1 


GeV- 2 


GeV" 4 


GeV" 6 


GeV" 8 


GeV- 10 


GeV- 12 


(0,7) 


+21257.2 


+146631 


+65353.0 


-26003.8 


+7140.37 


-997.823 


+56.4467 


(0,9) 


+10438.4 


+73176.5 


+32041.6 


-12776.8 


+3503.48 


-489.692 


+27.7109 


(0,10) 


+2821.90 


-122131 


-7771.38 


+ 178.016 


+20.9504 


-0.329242 


+0.0223758 


(7,7) 


-1326.63 


+57405.1 


+3778.94 


-61.4705 


-9.92262 


+0.105597 


-0.01443264 


(7,9) 


-1318.24 


+57047.6 


+3692.75 


-72.1696 


-9.82836 


+0.128021 


-0.0122690 


(9,9) 


-327.478 


+14173.1 


+901.859 


-20.6586 


-2.43126 


+0.0382081 


-0.00259668 



Table 14: Coefficients of the polynomial functions Ju j) entering F^. 



completeness), P5 and P$, while the different non-zero Ju t j) are either shown in table 
[TJ| or given by 



J(7',T) 


~ J (7,7), 


J(7',9') = —J{7,W) — 


— J(7',9) = 


^(7,9), 




J(9',9') 


= ^(10', 10') 


= ^(10,10) = <^(9,9), 


J(10,10') = 


: J(9,9'), 




J(0,T) 


= ~J(0,7), 


J{0,9') — —J(0,9)i 


J(0,W) = 


— ^(0,10), 




J(7,T) 


= — 2J(7,7) 


, and J( 9)9 /) = - 


"2 J(9,9), 




(2.45) 



rendering 20 entries Juj) different from zero entering F^. 

Therefore, the value of the integrated polarization fraction (Fl) in the SM can 
be computed theoretically using our inputs to get 

F[ M = 0.725i°;°37- (2-46) 
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2.6 B{B S -»■ 

The branching ratio of B° s — > n + fi~ in presence of only NP axial operators (relevant 
to this analysis) is given, at leading order, by [§, 



f 2 r> 2 / 4m 2 

B(B S /xV)Ui = m Bs r B jy t6 y t :| 2 m^ /l - -5^10 " CV| 2 (2-47) 

Using the inputs in table [I] and ^| we get the SM prediction 

B(B S -> ^ + /0 SM = (3-44 ± 0.32) • 1(T 9 , (2.48) 

which is one order of magnitude smaller than the most recent experimental averaged 
upper bound, obtained at the 90% confidence level in ref. [|3^] 12 : 



B(B S -> /i + ^) cxp < 3.2 • 10" 8 . (2.49) 

Eq. ( |2.47 ) can be used to compute a semi-numerical expression for this observable 
that will impose constrains in the (SCw, 8C\q>) plane (see Fig. |I]), 

B(B S -»■ = 1.8525 ■ 10~ 10 [| - 4.3085 + 5C W - 5C W \ 2 ± 1.7274] . (2.50) 

We have employed eq. Q2.50|) to check that the values of 8C\q and SCw used in 
Scenarios B and C (see below) were compatible with the constraints coming from 
B(B S — > fi + fi~). Since the experimental upper bound is still much larger than the 
SM prediction, no further cuts in the parameter space of Wilson coefficients have 
been found. 



3. Results 

In this section we obtain the allowed regions systematically from all the observables 
discussed previously. Since we aim first at illustrating how much our conclusions 
vary depending on the precise framework adopted to analyse the data, we will not 
adopt a sophisticated statistical approach (see refs. [[)7], [58| for examples of such 
approaches in similar contexts), and we will stick to a scanning approach, combining 
the 1 a theoretical and experimental ranges for each observable linearly to draw 
the corresponding constraint. For instance, if an observable Xj has been measured 
experimentally X t ± 8Xi and has the theoretical prediction Y^SCj) ± 5Yi, we draw 
the projection of the region corresponding to the constraint |JQ — Yi\< (5Xi + 5Yj). 

12 The LHCb Collaboration has just released a paper |33| where the upper limit on the branching 
ratio is set to B(B S — > < 5-6 • at 95% confidence level for an integrated luminosity of 

37 pb^ 1 . Since this upper bound is larger than the one obtained by the CDF collaboration [^2j we 
are not using it in this work. 
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Figure 1: Constraint imposed by B(B S — > to the values of the Wilson coeffi- 

cients in the (<5Ci , SC W ') plane. 

3.1 (C 7 ,C r ) plane 

As discussed in the introduction, we focus first on the C 7 , C 7 * plane, which will be the 
starting point of our discussion. Therefore, we consider the three Class-I observables 
which only depend on the electromagnetic operators C 7 ,C 7 ', leading to Fig. ^. If 
one considers only B{B — > X s ^) (ring in Fig. and S , ^» 7 (cross in Fig. |2|), four 
regions remain allowed: the SM one sitting around the origin, the "flipped- sign" 
solution discussed in the introduction around (5C 7 ,5C 7 /) = (0.9,0), and two non 
SM-like solutions with 5C 7 ~ 0.35 and 8C 7 i around ±0.5. The flipped-signed solution 
does not correspond exactly to C| — > —C^ E (and C 7 > ~ 0), due to interference terms 
between the electromagnetic operator and the four-quark operators in the observables 
considered here. The discriminating power of the isospin asymmetry in B — > K*j 
is quite obvious at this stage, as it discards this flipped-sign solution at 1 a without 
requiring further assumptions concerning NP for other operators. To recover this 
solution one needs to enlarge both theoretical and experimental uncertainties up 
to 1.59a. In our language, we disfavour this solution, working at la, on the sole 
basis of Class-I operators, contrary to ref. 0, [L8] which needed Class-Ill quantities 
[B(B — > X s £ + £~)} and thus obtained conclusions with more restrictive assumptions 
concerning the manifestations of NP. 

We will use the three identified black regions in Fig. ^] as the reference or primary 
regions: 

• the region around (5C 7 ,5C 7 /) = (0,0), referred to as the "Central" or SM-like 
solution; 
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Figure 2: Class I observables at 1 a: Aj (solid blue region with a white disk), B(B — > 
X s ^) (orange ring) and Sk* 1 (red cross). The three disconnected regions allowed by 
the intersection of these three observables are depicted in black. The SM value is 
given by the crossing of light gray lines at (SC7, SC?) = (0, 0) point. All plots of 
Wilson coefficients are taken at fib = 4.8 GeV. 

• the upper region around (<5C 7 ,<5CV) = (0.35,0.45), referred to as the "Upper" 
region; 

• the lower region around (SC7, 5Ct) = (0.30, —0.45), referred to as the "Lower" 
region. 

The last two regions will be commonly called non SM-like solutions in the following. 
These regions constitute the starting point to study the impact of Class-II and Class- 
Ill observables under the three different scenarios (A, B and C) presented in the 
introduction, each more general than the previous one. 

It is important to remark that the two non SM-like primary regions of Fig. |2] 
contain an interesting subset of solutions for CV with a flipped sign with respect to 
the SM. These solutions are characterised by a small modulus of CV and the addition 
of a larger contribution from C71 to get agreement with data. 
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Figure 3: Constraint from Class-Ill observables B(B — > X s fi + fi~) (left), Afb (mid- 
dle) and Fl at 1 a in the (£CV, £CV) plane in Scenario A together with the three 
(black) regions allowed by Class-I observables. The magenta circle centered at 
(0.25, —0.40) on the first plot indicates the tiny allowed region in this Scenario A. 



3.2 Scenario A 

Let us start with Scenario A. If we consider the Class-Ill observables B(B — > 
X s [i + fi~), A FB and F L for B — > K*fi + fi~ in the low-q 2 region, we obtain the con- 



straints shown in Figs. ga|, [3l] and [fc] respectively. One observes that the three ob- 
servables favour different regions of the (Cj,Cr) plane: the inclusive decay favours 
the SM region and a very small subregion inside one of the non-SM like solutions, 
whereas (as expected) the forward-backward asymmetry would favour the flipped 
sign-solution (had it not disappeared due to the isospin asymmetry) but also the two 
non-SM like solutions. The longitudinal polarisation would agree with all the regions 
(cutting only a very small part of the flipped-sign solution region). 

We see that Scenario A yields somewhat contradictory information from the var- 
ious observables concerning which region in the (C-j, Cf) plane should be preferred. 
There is actually only a very small region in perfect agreement with all the observ- 
ables measured (Class I and Class III), around 5C 7 ~ 0.25, 5Cr ~ —0.40 highlighted 



with a magenta circle in Fig. Ral corresponding to the intersection of the lower black 



region with the B — > X s fi + fi~ constraint. It makes therefore sense to extend the 
set of operators potentially affected by NP and to consider Scenario B, including 
also New Physics in Co, and Ciq. Before leaving Scenario A, it is very interesting to 
compute the values for the (Class-II) observable A}p that is not yet measured, and 
turn it into a prediction. Figure f| illustrates the prediction for this observable as a 
function of q 2 for the small set of points allowed by Scenario A. This leads to a very 
precise prediction for the variation of A^r (including error bars) with q 2 . Notice that 
all the curves included in this region exhibit a zero in a range between 1 to 1.6 GeV 2 
which is controlled, at LO, by the same equation that fixes the position of the zero of 
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Figure 4: Prediction for A T under Scenario A (lower pink band), including error 
bars for all the allowed New Physics curves. The band around zero corresponds to 
the SM prediction. 



A FB ||, |1. Given the small value of Cf(p b ) ~ -0.29 + 0.25 = -0.04, the position 



of this zero is shifted to the left with respect to the SM. Finally, another important 
prediction of this scenario is that A^\q 2 ) would clearly prefer negative values, due 

(2) 

to the negative value of Cy ~ —0.4. Therefore, a measured value for A\ different 
from the narrow prediction given here would be enough to rule out this scenario. On 
the contrary, a measurement consistent with this prediction would make Scenario A 
the most plausible one (compared to the other scenarios), and furthermore, would 
signal clearly the presence of right-handed currents in radiative decays. 

3.3 Scenario B 

In case of Scenario B, the regions permitted by the Class-Ill observables B(B — > 
X s fi + fi~), A FB and F L in Figs. ^a[ [TE|, [k] become extended to the whole plane, 
and thus are not constraining anymore either CV or Cy. In this scenario, the three 
primary (black) regions in Fig. |2] allowed by the Class-I observables are compatible 
with all the Class-Ill observables considered and become the allowed region for Cj 
and Cji in this scenario. This obviously does not mean that the observables of class 
III mentioned above do not provide any constraint on NP, just that these constraints 
are not visible in this particular subspace of NP parameters. As emphasized in 
the introduction, the (C 7 ,C^) plane is a summary that does not provide the full 
information on NP. It is thus interesting to turn to the (Cg, Cio) plane. Figs. [[a|, 



and |5c] are obtained taking the values of the (now) permitted three primary (black) 
regions in Fig. || for (CV, Cy) and determining the values of Cg and Cio that are then 
allowed for B(B — > X s /i + /i~), A FB and F L , respectively. 
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It is quite interesting to notice that the region excluded by Fl is very close to 
the central region excluded by B(B — > X s ^i + ^r). This is more striking once all 
constraints from the three observables are overlapped in one single Fig. |6|, where 
only two regions (shown in black) are allowed by all constraints. The nature of these 
two areas can be understood by in the following way: 

• SM region: the region centered at the origin corresponds to deviations from 
SM values for (Cg, Cio) keeping the same sign for these coefficients as in SM; 

• flipped- values region or non-SM region: this solution contains a subregion with 
opposite sign values for Cg and Cio with respect to the SM ones. 

The existence of these two regions can be understood from the fact that most 
of the observables have an approximate symmetry consisting in changing the sign 
of C 9 , C w altogether, as long as C 7 or C 7 > remain small (see, for instance, the large 
recoil expression for A FB in eq. ( |B.47| ) of Appendix |B.5.4| with Cg/ = C W i = 0). 



We checked that each of the three primary (black) regions in the (Ct,Cy') plane 
yield Class-Ill constraints in the (Cg, Cio) plane that cover the two regions in Fig. 6 
almost entirely. It implies that the two regions in (Cg, Cio) plane exist independently 
of the precise values for C7 and C7/, as long as any of the latter remain small and 
have a limited impact on the leptonic observables. In our framework, this smallness 
is indeed ensured by the constraints in (C7, C7/) plane coming from B(B — > X s ^y). 
It is interesting to provide predictions for the (still not measured) asymmetry 

(2) 

Aj, , using as inputs the WCs associated to the three black regions allowed in (C7, C7/) 
plane, together with the corresponding set of values in the (Cg, Cio) plane (two black 
regions). This is shown in Fig. |7|. We can see there that the large allowed areas for 
(Cg,Ci ) lead to wide bands in A^\q 2 ). The Upper non-SM like (C7, C7/) region 
associated to the SM-like (Cg, Cio) area gives a clear prediction for the sign of A 



(2) 
T ) 



which is just opposite to the one preferred by Scenario A. Also the Central (SM-like) 
(C7, C7/) region associated to the non SM-like (Cg, Cio) area (Fig. 7a) and the Lower 
(C7, C7/) region associated to the SM-like (Cg,Cio) area (Fig. 7f) yield constraints 

(2} — 

on , though less stringent than those in Fig. 7e. 

In conclusion, in this scenario the upper region of (C7, C7/) with the correspond- 

(2) 

ing SM-like region for (Cg, Cio) could be discriminated clearly only if the sign of A T 
would turn out to be negative, as predicted by Scenario A. Besides, high-g 2 measure- 
ments, not included in the present analysis, could shrink the allowed (Cg, Cio) region 

(2) 

and thus reduce the range of possibilities for A^' in this scenario. 
3.4 Scenario C 

Finally, we could imagine that the previous constraints did not overlap as nicely as in 
Fig. |6|. We would then turn to Scenario C, allowing for chirally-flipped semileptonic 
operators. For (8C7, 6Cr>), we take all the model-independent allowed values from the 
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Figure 5: Constraint from Class-Ill observables B(B — > X s fi + fi~) (left), Afb (mid- 
dle) and F L (right) at 1 a in the (6Cg, 5C W ) plane in Scenario B. The region shown 
is compatible with the constraints on 6C 7 and 5C 7 / imposed by Class-I observables. 



three regions of Fig. 0. Among all the constraints considered previously from Class- 
Ill observables, only BiB — > X s fi + fi~) still provides a constraint on the semileptonic 
(primed and unprimed) Wilson coefficients. Indeed, when NP contributions in Cg> 



and Cio' are also considered, the empty region in the middle of Fig. ^ gets filled 
up but the minimum and maximum values of 5Cg and 8C\g allowed do not change 
perceptibly, as can be seen in Fig. |8a]. In Fig. |E we show the allowed region in the 



(5Cg>, 5C W >) plane in the same scenario. It is not very surprising to obtain such oval 
shapes in the various planes of interest, since it corresponds to the projections of the 
quadratic (elliptic) constraint given by eq. ( p. 24 ). In conclusion in Scenario C, the 



allowed region for (6Cr,5Cr) is given by the three black regions in Fig. 0, and the 
corresponding ones for the planes (SCg, SC W ) and (SCg>, 5C W >) are given by Figs. |Sa 
and ^ respectively. 

(2) 

We have not given the predictions for under this scenario, as the extra 
freedom provided by Cgi and Cio' is likely to fill the whole parameter space available. 

3.5 2 a constraints 

When the uncertainty in both theoretical and experimental results is increased to 
2 a, the regions allowed in the (<5CV, 5Cr) plane are enlarged, as B(B — > X s j), Sx'-y 
and Aj yield larger overlapping regions. More importantly, the whole region corre- 
sponding to the "flipped-sign" solution is no longer excluded by Class I observables 



(see Fig. ^). We have followed the procedure explained before and used the re- 
sulting four disconnected regions to explore the behaviour of Class-II and Class-Ill 
observables under scenarios A, B and C. 

In Scenario A, B(B — > X s fi + fi~) excludes the whole "flipped-sign" solution re- 
gion, a sizeable portion of the upper region and small part of the lower one, as shown 
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Figure 6: Overlap of the constraints from Class-Ill observables B(B — > X s fi + fi~) 
(green ring), Afb (upper and lower "hyperbolic-like" brown regions; see Fig. 5b) and 
F L (dark gray area with a central inlet) at 1 a in the (6Cg, SCi ) plane in Scenario 
B. The constraints imposed by their intersection are shown as two black regions. 



in Fig. EJEL whereas neither v4fb nor -Pl provide further constraints, since they fill 



the whole of the (5C 7 , 5Cy) area explored. 

Next we move to Scenario B and include possible NP contributions to (SCg, SC W ), 
as depicted in Figs. pLQa| and |10b| . The region allowed by B(B — > X s fi + becomes 



enlarged by about a 40% with respect to the 1 a plot and the central region, pre- 
viously forbidden, becomes filled altogether. In this scenario, Afb does not provide 
extra constraints but maintains an excluded central zone, although much reduced 
in area. Fig. |TT| shows (in black) the regions allowed by the overlapping of these two 
observables. Moreover, the "flipped sign solution" for the (5Ct,SC 7 ') plane is now 
allowed under this scenario and the following one. 

We come finally to Scenario C. Besides (SCg, SCig), we must also allow for NP 
in the Wilson coefficients Cg> and Cio>, while (5CV, 5Ct) remain confined to the four 
black regions of Fig. |9a[ B(B — > X s fi + fi~) is again the only observable that imposes 
constraints in the Wilson coefficients related to Oi and Oi> (with i = 9, 10) as shown 
in Figs. [I2aj and |12b. 
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Figure 7: Prediction for A^\q 2 ) corresponding to the "flipped- values" region in 
(Cg, Cio) plane (first row of plots) and the SM-like region in (Cg, C\q) plane (second 
row) in Fig. 6. Each column corresponds to SM-like Central region (left), non SM- 
like Upper region (center), non SM-like Lower region (right) for the (C 7 , C 7 i) plane 
allowed regions in Fig. 2. 

3.6 Generalization to extended frameworks 

Let us assume, for instance, that we want also to include contributions from scalar 
operators (like those defined in ||). Consequently the scenarios will also be enlarged: 
Scenario A (0 7 , O r ), B (0 7 , O r , Cg, O w ), C (0 7 , O r , scalars), D (0 7 , 7 >, Cg, O w , O v , 
O w i), E (0 7 , 7 r, Og, Oio, scalars), F (all operators). We would then proceed again 
along the same steps as before, up to certain changes: 

1. We classify again the observables according to this new framework. This may 
move some observable from Class-I to higher classes, because they have sensi- 
tivity to scalars, like the K* / y observables Aj or Sk* t Only B(B — > X^) will 
remain. 

2. We determine the new reference region for C 7 and C 7 i defined by the (now 
reduced) set of Class-I observables. The new primary regions will be larger 
than in the previous framework because some observables are not included in 
the new Class-I. 

3. At this stage, and working in Scenario A, it is interesting to define two types 
of Class-II observables, Class-IIa, only sensitive to dipole, semileptonic and 
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Figure 8: Constraints from Class-Ill observable B(B — > X s fi + fi~) at 1 a in the 
(SC 9 , 5C W ) and (SCy, SC 10 ') planes in Scenario C. The regions shown are compatible 
with the constraints on <5CV and 5Cr imposed by Class-I observables. 



chirally flipped (our observables in Class II of the previous framework will be 
here) and Class-lib, only sensitive to dipole operator (and its chirally flipped 
counterpart) and scalars. These observables may shrink the new reference 
regions, leading to allowed regions of different shapes for Class-IIa and Class- 
lib. If we add now Class-Ill observables with sensitivity to the whole list of 
operators in the framework, this will generate a further cut on the primary 
region. If the same set of observables as in the previous framework has been 
included it is clear that, even if re-classified, the allowed region under Scenario 
A will be exactly the same as in the previous framework, even if the primary 
regions are different. 

4. The main differences arise when dealing with the rest of scenarios. We should 
repeat the same analysis under Scenario B till F. It is clear that Scenario B 
and Scenario C, for instance, may select different subregions inside the primary 
regions, and that Scenario D will enlarge the region for Scenario B, and the 
same will happen between scenarios E and C. Finally Scenario F will cover all 
previous ones, defining the largest allowed subregion inside the primary regions. 
This region maybe larger that in the previous framework (since more freedom 
in the value of the WC already studied is provided by the introduction of scalar 
contributions). 

5. This systematic procedure that subdivides the primary regions in different sub- 
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Figure 9: On the left, 2 a constraint from Class-I observables: Aj (background solid 
blue region with two white disks -partially hidden-), B(B — > X s '~f) (orange ring) and 
Sx'-y (red cross). The three disconnected regions allowed by the intersection of these 
three observables are depicted in black. On the right, 2 a constraint from Class-Ill 
observable B — > X s/ u + /i~. The SM value is given by the crossing of light gray lines 
at (5C 7 , 5C V ) = (0, 0) point. 

regions may help to disentangle the importance of each set of operators: dipole, 
semileptonic, chirally flipped, scalar, when confronting theory with data. In 

(2) 

particular, certain observables like and its generalization, may discriminate 
between the different subregions. 

This procedure can be generalized to other frameworks following the same steps. 
Defining intermediate steps between the dipole-only case and the full-fledged scenario 
for the introduction of New Physics helps in understanding the importance of the NP 
contribution for each observable. In the present paper, we have restricted ourselves to 
the framework where NP arises in dipole, semileptonic operators and their chirally- 
flipped partners. 

4. Discussion and outlook 

We have exploited the (SC 7 , 5C 7 i) plane as a starting point to investigate the pattern 
of NP in the Wilson coefficients for radiative AB = 1 transitions. We have defined 
several classes of observables to help us in this task, selecting only observables with a 
good theoretical control over hadronic uncertainties (or a significant discriminating 
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Figure 10: Constraints from Class-Ill observables B(B — > X s fi + fi~) (left) and Fl 
(right) at 2 a in the (5Cg,5Cio) plane in Scenario B. 



power for our NP scenarios) and providing numerical expressions for these quanti- 
ties as functions of ^CVj'^g^io.io'- We defined reference regions for (5C 7 ,5C 7 r) from 
Class-I observables, then studied several scenarios of NP involving chirality-flipped 
operators with the help of Class-II and Class-Ill observables. 

As far as the theory and experimental errors of the measured observables remain 
inside the 1 a range we can draw the following conclusions. Scenario A, where only 
(C 7 ,C 7 >) receive large NP contributions, is a predictive scenario. Class-I observables 
provide three different regions (one corresponding to the SM case, two other ones 
with almost vanishing C 7 values and large C 7 > value). Once Class-Ill observables are 
included only a very small subregion (inside one of the two non-SM like regions) is 
allowed if we keep all the constraints at 1 a . Consequently, only those theories that 
can provide values for (CV, G 7 i) ~ (C 7 + 0.25,-0.4) are compatible (within la) 
with current (Class-I and Class-Ill) data, due to the interplay between the inclusive 
decay B — > X s fi + fi~ and the forward-backward asymmetry Afb- Notice that the 
SM is not one of such theories. This motivated us to enlarge the set of operators 
where NP contributions can be sizeable, leading to constraints on the semileptonic 
operators. Scenario B constitutes the first extension, allowing for NP in CV,7',9,io- 
In this case, the previous constraints from Class-Ill observables are transferred from 
the (SC 7 ,SC 7 ') plane to the (5C 9 ,5Ci ) one. There are two distinctive regions al- 
lowed, corresponding to the SM solution, but also to a flipped-value configuration, 
where Cq and C%q have some values opposite to the SM. It is interesting to notice 
that B — > X s fi + fi~ and Fl exclude almost the same central area in the (SCg, SCiq) 
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Figure 11: Overlap of the constraints from Class-Ill observables B(B — > X s fi + fi~) 
and Fl at 2<T in the (5Cg,SCio) plane in Scenario B. The constraints imposed by 
their intersection are shown as a black region. 

plane. Scenario C (with NP in CV,7',9,9',io,io') would be an interesting extension if 
the previous experimental constraints shift in the future, or if the measurement of 
the (Class-II) asymmetry Arj? shows a discrepancy with the pattern of Wilson coef- 
ficients exhibited in Scenario B, once more data and constraints have been added. 
The (Class-I) constraints on (CV, CV) remain unchanged with respect to Scenario B, 
whereas Class-Ill observables provide only limited constraints on the largest set of 
Wilson coefficients considered. Currently, only B — > X s fi + pT provides constraints 

On 

We have also indicated how the (Class-II) asymmetry A T ' gives a very precise 
prediction for Scenario A, that can be used either to confirm it or to rule it out. 
It may also help, depending on its sign, to discriminate among the allowed regions 

(2) 

in Scenario B. exhibits a strong sensitivity to the allowed regions for (Cg,Cio); 
further cuts in these regions using high-g 2 measurements, will improve the predictive 

(2) 

power of A^' in this scenario. Under Scenario C, there is too much freedom with all 
WC switched on to be able to cut on precise regions as it happens for most of the 
other observables. 

We also have shown that Class-I observables alone allow us to dismiss the flipped- 
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Figure 12: Constraints from Class-Ill observable B(B — > X s fi + fi~) at 2 a in the 
(SC 9 , 5C W ) and (5C 9 ', SCw) planes in Scenario C. The regions shown are compatible 
with the constraints on 5CV and 5Cr imposed by Class-I observables. 



sign solution at 1.59a, even in a NP scenario much more general than in ref. [[/J, 
allowing for NP in dipole and semileptonic operators, but also in their chirally-flipped 
counterparts. We achieved this by trading the Class-Ill observable B — > X s fi + fi~ 
(considered in ref. J7), and sensitive to many NP contributions apart from those in 
the dipole ones) for the Class-I isospin asymmetry in B — > K*^ (even though the 
theoretical control on hadronic uncertainties is less satisfying for this observable). 
A summary of the maximal and minimal values of the WC analyzed in the 



different scenarios is provided in table [15 



13 



In ref. [|ll|, an analysis of various NP contributions was considered, allowing 
either for New Physics in (CV,CV') (both of them being real), or C\o (considered as 
potentially complex). In particular, our findings concerning Scenario A (NP only 
in C-j and C?>) are in agreement with Fig. 2 in ref. ]TB| concerning S#* 7 , as well 
as the fact that the flipped-sign solution is excluded (even though the conclusion 
is based on different observables). However, the other scenarios discussed in [|18 



considered NP entering in one Wilson coefficient at a time, and thus provide only 
a particular section of the parameter space of Wilson coefficients. Another related 
study was performed in ref. [|5J|, where B — > K*£ + £~ at large and low recoil (which 
was not considered here) was combined with B — > X s £ + £~ to study the (Cg, C w ) 
plane, considering CV = ±Cj M . This led to two regions in (Cg, C*io) similar to the 



13 For the internal 4-d and 6-d correlations involving 4 WCs (Scenario B) and 6 WCs (Scenario 
C) we can provide a datafile with the correlated points upon request. 
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SC r (p b ) 


SC 9 (/j b ) 




$C 9 > (fi b ) 


5C 10 > (fi b ) 


Overlap of the 1 a constraints 


Sc. A 


[0.244, 0.274J 


[-0.417, -0.39J 














Sc. B 


[0.346,0.385] 
[-0.056,0.016] 
[0.235,0.385] 


[0.435,0.501] 
[-0.114,0.027] 
[-0.489,-0.39] 


[-9.75,-0.5] 
[-3.75,3.5] 


[4.75,10.5] 
[-1.75,3.5] 








Sc. C 


[0.346,0.385] 
[-0.056,0.016] 
[0.235,0.385] 


[0.435,0.501] 
[-0.114,0.027] 
[-0.489,-0.39] 


[-10,3.5] 


[-1.5,10.5] 


[-8,8] 


[-6,6] 


Overlap of the 2 a constraints 


Sc. A 


[0.262,0.586] 
[-0.083,0.076] 
[0.124,0.475] 


[0.381,0.531] 
[-0.225,0.105] 
[-0.519,-0.306] 














Sc. B 


[0.262,0.646] 
[-0.083,0.076] 
[0.775,0.97] 
[0.124,0.481] 


[0.381,0.534] 
[-0.225,0.105] 

[-0.12,0.3] 
[-0.519,-0.306] 


[-13.4,4.5] 


[-2.5,11.4] 








Sc. C 


[0.262,0.646] 
[-0.083,0.076] 
[0.775,0.97] 
[0.124,0.481] 


[0.381,0.534] 
[-0.225,0.105] 

[-0.12,0.3] 
[-0.519,-0.306] 


[-13.5,4.6] 


[-2.6,11.5] 


[-9,9] 


[-7,7] 



Table 15: Summary table of the maximum and minimum Wilson coefficients values 
allowed by the three different scenarios within our framework. The table is organized 
in three independent blocks corresponding to the pairs (5C 7 ,5C 7/ ), (SC 9 ,5C 9 i) and 
(5Cio,5Cio>) respectively Notice that the correlations between different WCs are 
more complex than those summarised in this table. In order to recover the exact 
2d-correlations, one should look at Fig. 3a (Scenario A), Figs. 1, 6 (Scenario B) and 
Figs. 1, 6, 8a, 8b (Scenario C) at la, and at Fig. 9b (Scenario A), Figs. 9a, 11 
(Scenario B) and Figs. 9a, 11, 12a, 12b (Scenario C) at 2 a. 

ones obtained in our case, however smaller partly due to the additional constraints 
put on C-j (and Cji) in this reference. 

In ref. ||57|| , a global analysis of AB = 1 observables was performed in a minimal 
flavour violating framework that included the possibility of sizable scalar contribu- 
tions (but no chirally flipped operators). The combination of the various observables 
was performed using a Bayesian statistical approach. Even though the inputs and the 
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underlying assumptions concerning the structure of NP are different (scalar versus 
chirality- flipped operators), we observe some common features. Two different regions 
for (C 7 , Cg, Cio) are allowed, corresponding approximately to a change of sign for the 
Wilson coefficients (Fig. 1 in ref. [0). Once NP is allowed for (C 9 ,Ci ) (Scenario 
B), there is a ring-like constraint from BiB — > X s £ + £~) in the (Cg, Cio) plane, with 
only two regions surviving once the forward-backward asymmetry Aps is included 
(Fig. 4 in ref. [J57| ). This is in basic agreement with our own plots, even though 
we should highlight that the non-SM region in the (Cg, Cio) plane corresponds to 
different allowed values for the electromagnetic operators: in ref. |57|, this region 
corresponds to the SM and the "flipped-sign" solution (C7 ~ — Cf M , Cy ~ 0) dis- 
favoured by B — > X s fi + fi~ in their framework, whereas our region corresponds to the 
SM solution and to the flipped- value regions where C 7 ~ and |C 7 /| ~ |C 7 |. 

Our approach could be extended to other, more involved, scenarios of New 
Physics, including contributions to the chromomagnetic, scalar and/or tensors op- 
erators as explained in detail in Sec |3.6| , allowing us to assess the impact of each 
observable in a controlled way. Such a task is left for future work. 
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A. Inputs 



We have followed the discussion in refs. 



E31 m 



concerning the matching 
and the running of the Wilson coefficients from the high scale fio = 2Mw down to 
the low scale fib = 4.8 GeV. We were able to reproduce at the 1% level the tables 

[J5j for 



3, 4 and 5 in ref. 



21 (apart from . 



(22) n (22), 



10 



and the table 5 in ref. 



the Wilson coefficients, providing a check that we control the scale dependence of 
the Wilson coefficients accurately. Contrary to other analyses in the literature, we 
have expressed the deviations from the SM Wilson coefficients at the low scale fit, 
around 4.8 GeV. However, the evolution from fi Q to fit, can be determined as the 
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linear combinations: 



5C 7 {fi b ) = 0.575 x 6C 7 {fio) , 
5C 9 (/i b ) = 1.021 x 5C 9 (hq) + 0.008 x 5C 10 (fx ) 
6C w {fi b ) = 0.008 x «yC 9 (/io) + 1-038 x 5C 10 (n ) 

Several schemes have been used to define the quark masses: 



(A.l) 



For m t and m c , we used the MS scheme scheme at the required scale (respec- 
tively /j,o and m c 
in refs. [26, 30]. 



We convert m? t ole into m^ s using the conversion formulae 



For mb, two different masses are needed: the mass in the IS scheme (or an 
equivalent scheme with infrared subtraction) is required whenever the 6-quark 
is close to the mass shell, whereas the pole mass is used for normalisation 
purposes as well as for loop computations where the 6-quark is off-shell. Foil- 
wing ref. [22, [35|, we take the value of m\ s obtained from fits to hadronic 
and leptonic moments of the differential branching ratio for the inclusive decay 
B — > X c iv |[T) , and we determine the pole mass using the conversion formulae 



in ref. 521. 



For m s , we use the strange quark mass in the MS scheme, taken at the scale 
fib- We are aware that there is an ambiguity in the scheme and scale chosen 
for this mass (this ambiguity would be resolved by going to higher orders in 
perturbation theory, which are not included in the present analysis). We used 
m s /mb both in the MS scheme to evaluate the SM value of C 7 > (however we kept 
the ml° le normalisation to determine rh s = m s /mb needed for B — > X s £ + £~). 



The running of the quark masses in the MS is performed following ref. |22fl . The 
strong and electromagnetic coupling constants are determined by their value at Mz, 



and their running is given by the equations in ref. [22 . 



B. Extension to chirally-flipped operators 

B.l B -»■ X sl 

The branching ratio for B — > X s ^ for a photon energy larger than Eq =1.6 GeV can 



be written as [35 



B{B -> X sl ) E ^ >Eo>SM = B(B -> X c ez7) 



V, 



cb 



Qa cn 
Ctt 



■[P(E ) + N(E )), (B.r 



- 38 - 



where 



C 



ub 



cb 



T{B ->■ X c eu) 



P(Eo)= Cf^Ct^K^Eo^). 

ij=1...8 



(B.2) 
(B.3) 



Concerning B — > X s 7, we were able to reproduce, not only the central value and 
uncertainty for the branching ratio, but also the results from the three different in- 



terpolation procedures and the scale dependence on /x and fib described in ref. |35 



as well as the dependence on CV j8 at the scale /i in eq. (29) of ref. [Q. The con- 
tribution from the chirally-flipped operator should have the same structure as 
the SM operator 7 and there are no interferences between the two contributions, 
leading to an additional contribution to eq. ( |B.1| ) of the form: 



P(E ) P(E ) + (C r ) 2 [l + a a (p>)K$ + K. 



'(2)1 



77 J ' 



(B.4) 



where are the coefficients of the perturbative expansion of the kernel K 77 (E , jj). 
B.2 £> — )• K*j isospin asymmetry 

Concerning the isospin asymmetry, we reproduced the central value of the isospin 
asymmetry quoted in ref. [[31] , following the formalism discussed in ref. |39 1 : 



Aj[B ->■ iT 7 ] 



SM 



Re[tf(0) - #(0)] 



SM ■ 



J q,SM 



(0) 



12n 2 f B e q 
m b C 7 £ ± (0) 



m B 



Ki{0) + 



fK* m K* 

6\ B m B 



i^(0) 



(B.5) 
(B.6) 



where C 7 = C 7 S + 0(a s ) includes NLO corrections to the amplitude for B — > K*j, 
computed in ref. [Q. In K^ify, we have included the Cabibbo-suppressed power 
corrections discussed in App. A. 3 in ref. and neglected in ref. |39j, performing 
the replacements 



.J_(c) 



4% 



^,™2 



um 



BJ 



C 2 



Xt 
Ci 
6 



[h(s,m c ) - h(s,0)] 



(B.7) 
(B.8) 



following the notation in ref. 39 



Unfortunately, the hard-spectator scattering involving the chromomagnetic op- 
erator Og exhibits an endpoint divergence indicating a breakdown of QCD factori- 
sation. We follow refs. JH, 



to regularise the divergent integral 

•1— Afc/m-B 



du — Y (1 + pe 



i<t>\ 



du , 



(B.9) 
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Figure 13: Annihilation topologies involving operators C?i_6(left). Hard spectator 
interaction involving operator 8 (center) and Oi-q (right). 



where p is assumed to be smaller than 1 for our numerical estimations, and the phase 
is arbitrary. 

Once we add chirality-flipped operators, Or will contribute to the branching 
ratio of B — > K*j. It is not difficult to check that its contribution is the same 
as the one from Of, and that there are no interferences between the two contribu- 
tions. We will neglect the contributions from SM operators to the amplitude for a 
photon of right-handed helicity. On the other hand, the flipped operators consid- 
ered in the present paper do not contribute to the spectator interactions responsible 
for the isospin asymmetry (which are induced by the four-quark operators and the 
chromomagnetic operators). Therefore, the only change induced by chirality-flipped 
operators corresponds to modifying the normalisation, i.e., the denominator in the 
expression of the isospin asymmetry (at first order in isospin breaking) 



A A B -> K-;] = R f + %-^r ■ (B-10) 



B.3 S K * 



We define the decay amplitudes of B d mesons into K* and 'Jl(b) as i n IHI : 

A m = A(B° d K*° lL[R) ), A L{R) = A(B° d K*° lL{R) ). (B.ll) 



With the assumptions explained under eq. (p.llf ) and using eqs. ( |B.11|) , the mixing 



induced CP-asymmetry (S) and the direct CP asymmetry (C) can be written as 



9= 21m [r d (A* L A L + A* R A R )] 

' |^Il| 2 + |^| 2 + |A| 2 + |^| 2 ' 



\A L 


2 + 


A R 


2 - 


A L 


2 -\A r 


2 


\A L 


2 + 


Ar 


2 + 


A L 


2 + \Ar 


2 



(B.12) 
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where = e and <pd is the 5° — 5° mixing angle. 

In "naive" factorisation, the decay amplitudes of eqs. ( |B.11|) are given by 



A L 



4G F 

"71" 

AG F 

' 



(K* lL \<D$\B), 



\(s)M u) + \ (S 



W (Mt) 



(B.13a) 
(B.13b) 



and 



At 



Ar 



AG F 

"71" 

AG F 

"71" 



(Ai s) )*^SM + (^)* (4'!sm + C<? ) (K*1l\ {0?Y\B), (B.14a) 



»(*) 



(B.14b) 



where, we have used the short-hand notation introduced in eq. 

Mi) - Tlr® 

W'.SM ~~ „ W,SM 
iil'b 



(B.15) 



with q = u,t. We have taken the notation and definitions from ref. H4~0| : Cf 1 are 
coefficients, defined as a ratio of full form factors and soft form factors, that can be 
computed in QCDF (Cy is equivalent to C| ff at LO in a s whereas Cj vanishes). 
Setting Cfft SM = and taking real Wilson coefficients Cj S and Cy, the mixing- 



induced CP-asymmetry yields the simple tree-level expression in ref. |18| , |41| , |42 

sin (2/3 - arg (cf (0) C r ) ) . (B. 16) 



q(lo) 



-2\C r /Cf (0) | 

1 + |c 7 vcf (0) | 2 



Eq. (|B.16| ) determines the cross-shaped plot of Sk*~, in the (SC 7 , 5C 7 i) plane (see 
Figure |2|) to a very good degree of approximation. We checked that S^?^ allows us 
to recover, at 2 a, the shape of Figure 2 (left) in [|18[] using their input parameters. 
Notice, however, that our actual computation, used for the plots in the present 
article, is performed including NLO QCDF corrections. 

Some comments are in order here. On the one hand, the operators O^j are 
given by 

°7 {R) = T^rribsa^^^bF^, (B.17) 



16tt 2 " 2 

and generate the left- (right-) handed photons in the b — > 37 decay. Following 
refs. f43| , f44|| we express the matrix elements in eqs. (p,13b| ) and ( p,14b| ) in terms of 
the form factor Tf~* K *(q 2 ) as 

(K*(p,vhL( R) (q,e)\Oi {R) \B) = 

m h T^ K " (0) {e^e; V tp p q a ± i [(eV)(w) - (e*p)fo*g)]} 



57T 

e 



m b Tf^(0)S 



m b T^ A (0)S L(R) 



(B.18) 
(B.19) 
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Figure 14: Form factor Tf~* K *(q 2 ) in the — 6 GeV 2 energy range. The orange bands 
represent the full form factor with its associated errors given by the parametrisation 
in Appendix B.4 of |)4j] and the gray dashed lines depict (q 2 ) computed from 



£,±(q 2 ) using the large- recoil expressions in refs. [pi 

where Sl,r are the helicity amplitudes corresponding, respectively, to left- and right- 
handed photons and e^rj^) is the polarisation four- vector of the photon (K*). 

On the other hand, since the photon emitted in the decay b — > s'j is real, only 
the operators Oi,...^ of the weak effective Hamiltonian contribute to this process. In 
particular, those that build up C% s (see eq. (|2.5|) ) appear at while the rest 

of the operators enter the NLO QCDF corrections. Even though there is just one 
form factor contributing to this process, we have used the corresponding soft form 
factor computed by means of eq. ( p.41|) to be consistent with the QCDF formalism 
applied to both Aj and B^ — > K*°£ + £~ observables. This amounts to replacing 
Tf^ K *(0) — > £j_(0) in eqs. ( |B.19| ), which is indeed a very good approximation, as we 
can see in Fig. |1"4] . 

Contrary to ref. [40 we have chosen to keep the CKM-suppressed terms pro- 



portional to Xu ■ In "naive" factorisation, both and Cy SM vanish at LO in 
a s . If NP is absent, C$ vanishes, as we have split the df) SM helicity-suppressed 
7# terms already present in the SM (see eq. (|2.4j) ) from the Of NP contribution. 
Therefore, including NP in the decay amplitudes can be obtained upon the following 
replacements in eqs. flB.13b| , |B.14b|) : 



-)• 6C r , Cf -> 4*sm + 6C ^ ( B - 20 ) 
where, as said in the previous section, df = Cf [q) + 0(a 8 ) includes the NLO 



corrections to the decay amplitude B — > K*^ [p3| . Therefore, the replacement 
^,±(0)C^ — > T^ 1 j21, ^0|, |63| in the expressions above will be enough to account 



for these corrections in QCDF. Using this framework, we have computed the 0(a s ) 
factorisable and non-factorisable corrections to hard-spectator scattering diagrams, 
as well as to those diagrams that involve a B — > K* form factor We have 
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also included the power-suppressed weak annihilation and hard-spectator scattering 
contributions following |^0|, |S5[| ; the latter suffer from the same kind of endpoint 
divergence that we find in Aj, and they have been regularised by means of eq. ( |B.9| ). 



B.4B4 X s £+£~ 

The branching ratio for B — > X s £ + £~, normalised by B — > X c lv and integrated 
between 1 and 6 GeV 2 can be written in the following manner: 



dB{B X 



ISM 



ds 



B(B X c eu) 



v t * s v tb 



cb 



4 $u(s) 
C K 



m 



b,pole 



(B.21) 



where 



i<3 



Cf(fi)Cf*(ti) Yl M t A Mf*S AB + AH i: 



\A,B=7,9,10 



(B.22) 



and C has already been defined in eq. 

We were able to reproduce the central value and uncertainty of B — > X s £ + £~ , 
but also the dependence on 6*7^,9,10 at the scale /io in eq. (12) of ref. J22| (apart 
from the linear term in CV(/io) which is very sensitive to small changes in the input 
parameters) . 



We have modified the building blocks S following ref. |66j to include m s correc- 
tions and contributions from chirality-flipped operators in the following way 14 . 



For the functions involving only A,B = 7, 9, 10, we modified the functions to 
include m s -suppressed contributions to the phase space and to 0(a®) part. 

For the functions involving only A,B = 7', 9', 10', we took the same expression 
as their unprimed counterparts, profitting from the fact that the expressions 
are symmetric with respect to the change 75 — > —75. 

For the functions involving both a SM operator and a chirally-flipped one, we 
took the expressions from ref. [p6 , which include only 0(oP s ) contributions 



(contrary to the other functions that include also 0(a s ) and 0(l/m 2 ) correc- 
tions). 



14 We checked and agreed with the expressions in ref. j6f|, taking into account the fact that this 
reference uses a different definition of Oj and 0' 7 which mixes different chiralities, contrary to ours. 
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5*77 — S-ji-ji 
S79 = SV'9' 
S99 — Sg'g' 



AM 1 + 
N (l + 

N 



2m\ 
s 

2m 2 



-4s - 4(1 + ml) + 8(1 ^ + 0(a s , l/m 2 b ) 



■ 12[l-m 2 s -s + 0(a s ,l/m 2 b )} 



1 + 2m] - 2m\ + 2m\m\ + m\ 



(B.23) 
(B.24) 

(B.25) 



2m 2 (l -m 2 ) 2 
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-12m s )(l - m 2 + s) 



S 99 , = jV(-12m s )(s + 2m 2 
Sioio' = iV(-12m g )(s - 6m 2 ) 

with the phase space factor 



N = y/1 + s 2 + m 4 s - 2s - 2m 2 s ~ 2sm 2 s \ 1 - 



Ami 



(B.28) 

(B.29) 
(B.30) 



(B.31} 



For the quantities related to matrix elements Mf, we have taken the expressions of 
ref. |22| for the unprimed operators. The situation is much simpler for chirally-flipped 
operators since only three of them are to be considered: 

Mj' = a s K5 i>r , M? = (1 + a s K^ en {s))5^, Mf = 5 i>w . (B.32) 

The uncertainty attached to the central value in table || includes not only the un- 
certainties from the variation of the difference input parameters, but also a 5% error 



estimated in ref. |22| as the uncertainty from non-perturbative l/m fe -suppressed con- 
tributions. 

B.5 5 4 K* £ + £- observables 
B.5.1 General considerations 

The differential decay amplitude of the exclusive process Bd — > K*°£ + £~, with K*° — > 
K~n + on the mass shell, can be characterised completely in terms of the dilepton 
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pair invariant mass q 2 , which is embedded in the so-called angular coefficients, and 
the three independent angles 9i, 6k and (see Section 2.1 of ||). These angular 
coefficients J, are observable quantities that depend on kinematical parameters, real 
combinations of the six complex K*° spin amplitudes and the seventh transverse 
amplitude A t (in the presence of scalars an extra amplitude is required |30j ) . 

Within our framework, the spin amplitudes can be expressed in terms of the 
seven B — > K* form factors and the Wilson coefficients Ci of the weak effective 
Hamiltonian, that account for the short-distance interactions. Neglecting 0(a s ) 
corrections and using the effective Wilson coefficient associated to Oj (which includes 
the contributions from the four-quark operators Oi,„g), as well as the numerically 
relevant coefficients C 9 and C 10 associated to 9 and O w respectively, we find 0: 



A 



L,R 



NV2X 1 ' 2 
2mi> 



{(C 9 + C 9 ,)t(C 10 + C w )} 



iris + rriK* 



+ 



+ 



{Cf + Cf^q 



(B.33) 
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K" 



{(C 9 -C 9 ,)t(C w -C w )} 



A,{q 2 



m B - m K * 



+ 



2m h 



{Cf-Cf)T 2 {q 



(B.34) 



A L,R 



A t 
where 
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2m 



K* 



{(C 9 -C 9 ,)T(Cia-C 10 ,)} 



[m 2 B - m 2 K * - q 2 ){m B + m K *)Ai(q 2 ) - 



+2m b (C 7 
NX 1 ' 2 



XA 2 (q 2 ) 
Tn B + m K * 

° S Cf^Uml + Zm^-q^Uq 2 ) 
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(B.35) 
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q- 



(B.36) 

(B.37) 
(B.38) 

(B.39) 



We have introduced the Wilson coefficients corresponding to the chirally flipped 
operators O-?/, 9 i and Ow, so we consider only NP contributions stemming from 
the SM-like operators and their chirally-flipped partners (i.e. we assume there are 
neither scalar /pseudoscalar nor tensor/pseudotensor operators at work). 
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B.5.2 Soft form factors 



Concerning the B — > K* form factors, there are seven a priori independent hadronic 
form factors, encoding the non-perturbative long-distance interactions, that enter the 
B — > K* matrix elements, namely the vector current form factor V(q 2 ), the three 
axial current form factors A (q 2 ), Ai(q 2 ), A 2 (q 2 ), the tensor form factor Ti(q 2 ) and 
the pseudo-tensor form factors T 2 (q 2 ) and T 3 (q 2 ) [pf] . Although there are several 
computations of these form factors in the literature (see for instance ref. p7|), we 
have chosen the parametrisation in Appendix B.4 of ref. p3 to remain more conser- 
vative in the estimation of the uncertainties associated to the fitting coefficients that 
describe them. In the limit where the decaying hadron is heavy (as in Bd) and the 
recoiling meson acquires a large energy (Ek*), the form factors can be expanded in 
the small ratios Aq C d/^6 and Aq CD /E k *. Neglecting corrections of order Aqcd/^ 
and a s , the seven B — > K* form factors reduce to just two universal "soft" form 
factors and £y J|8], |68 . 



In this limit the K* spin amplitudes and A t acquire very simple forms which 
prove to be most useful to explain the symmetries between the fitting coefficients 
(F, G, H, I, J and K) of B d K*°£ + £- observables given in Sees. O and O [plj 



A L ^ R = V2Nm B (l - s) 



2m h 



(C 9 + C v ) =F (C 10 + c 10 o + ^(Cf + Cf) 



An' R = -V2Nm B (l - s) 



Irhh 



(Cg — C 9 >) =f (Cio — CioO H r^-(Cj G — Cy, tt ) 
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NrriB 



(C 9 - C 9 =f (C 
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c 



10' 



2m b (C 



'Cff 



'off 



Gel 
7' 



(B.40a) 

Z±(Ek*), 
(B.40b) 

Z\\(Ek*), 
(B.40c) 

(B.40d) 



Cio — Cio' i\\(E K * 

mK*Vs 

with s = q 2 /m 2 B and rhi = rrii/mB- 

The QCDF framework allows us to calculate the a s corrections to form factors 
and decay amplitudes up to the NLO PD), [ISl in a systematic way but, since 
we have no means of computing the 1 /mb-suppressed corrections, we decided to 
estimate them consistently using an ensemble method for the K*° spin amplitudes; 
an exhaustive discussion of all these issues can be found in sections 2.2 and 2.3 of 
ref. i. 

However, since QCDF uses only soft form factors and not the full form factors, 
we are restricted to the kinematic region in which Ek* ~ rrib (or equivalently, q 2 
tub)- Moreover, the longitudinal spin amplitude displays a logarithmic divergence as 
q 2 0, which signals the breakdown of QCDF for energies below 1 GeV 2 . Further 
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cuts are provided by the light (below 1 GeV 2 ) and J/ip (over 6 GeV 2 ) resonances. 
Thus, we have confined the analysis of A\ , Afb and to the dilepton mass range, 
lGeV 2 ^ q 2 < 6 GeV 2 . 



We obtain the soft form factors demanded by the QCDF framework [63], |4(| from 
the full form factors V(q 2 ), A x (q 2 ) and A 2 (q 2 ) using |, [18], ||, |K§ 



m B 



m B + m K - 
m B + m K 



2E 



K* 



■V(q 2 ), 

-A(g 2 ) 



m B - m K , 
m B 



(B.41) 



Our choice of ref. |53| with sizeable error bars compared to other possible determina- 
tions is guided by our aim to be conservative in our estimation of errors. Eq. ( |B.41|) , 
in particular, defines the value of the soft form factors at q 2 = from the values of 
the full form factors taken from ref. |M . 



B.5.3 The differential decay distribution and uniangular projections 

The angular dependence of the B d — > K*°£ + £~ differential decay distribution can be 
integrated out yielding, in terms of the K*° spin amplitudes, 



dT 
dq 



, l [(Z + Pl)(\M 2 + \A\\\ 2 + \A,\ 2 )\ + 
3m 2 

+ —Ji {\A t \ 2 + 2[Re(A ±L A* ±R ) + Re(A\\ L A* R ) + Re(A 0L A; R )} } (B.42) 



where we have defined A { Aj = A iL A* L + A iR A* R , with i, j = 0, _L, ||. 

The large uncertainties coming from the B — > K* form factors turn dT /dq 2 into 
a theoretically ill-controlled observable (as can be seen in Fig. |15D. However, since it 
appears only in the denominator of Afb and F^, and the corresponding numerators 
display the same kind of uncertainties correlated to those in dT/dq 2 , Aps and F^ 
become much better behaved observables (see Figs. |16a| and |16b| , and refs. 
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an in-depth discussion of this issue). A^\ on the contrary, is essentially free from 
this problem. 

As shown in refs. ||, |47]], a full angular fit can be performed on Bd — > K*°£ + £ + 
observables, but this will probably require more integrated luminosity than the one 
delivered by the end of the fist run of LHC [H7L B9|, [70]]. However, we can also 



integrate out two of the three angles of the K* differential decay distribution to get 
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3 4 



Figure 15: SM prediction for the differential decay distribution of B& — > K*°£ + £~ 
in the 1 — 6 GeV 2 energy range. The black line corresponds to the central value of 
dT/dq 2 . The wide gray band corresponds to the uncertainties associated to B — > K* 
form factors (according to the parametrisation in Appendix 4 of ref. @). Hadronic 
(orange) and AqcdA 7 ^ (green) uncertainty bands are barely visible. The central 
value compares well with Fig. 2 in ref. (note that the CP-averaged differential 
decay distribution d(T + T)/dq 2 is plotted there, so their central value curve is twice 
ours) . 




1 2 3 4 5 6 1 2 3 4 5 6 



q^GeV 2 ) q^GeV 2 ) 
(a) (b) 

Figure 16: SM prediction for Ape (left) and (right) in the 1 — 6 GeV 2 energy 
range. The color scheme used for uncertainties is the same as in Fig. 15. 



three single-angle distributions, which, in the massless case read 

P d0 = 27 ( 1+ 2( 1 " FL ) A T 2) cos20 + A im sin20j , (B.43a) 
(^Lsin 2 ^ + ^(l-F L )(l + cos 2 ^) + ^F B cos^ sin0,, (B.43b) 
-sin^ (2F L cos 2 ^ + (1 -F L )sin 2 ^) , (B.43c) 



V d9 K 4 



15 Massive terms are suppressed by m 2 ~ 0.011 GeV 2 so that their impact in absence of possible 
large scalar/pseudoscalar or tensor/pseudotensor NP operators is negligible. 
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where 

V =—2 and A ™ = dT — ■ (B.44) 

dq 



dq 



2 



(2) 

Since Afb, F^ and A T appear in all the expressions above, experimental data can 
be binned in q 2 and the corresponding fits performed on these bins. The value 
extracted from these fits is then a ^--weighted average of each parameter. Some 
strategies have already been devised to perform the binning in a way that allows to 



increase the statistics signal for some chosen observables |54], [7(J |7T . 

B.5.4 B d — > K*°£ + £~ observables at leading order in the large-recoil limit 

This section is devoted to the analysis of the relations existing between the different 
fitting functions, namely eqs. ( [T2CD , (gTglj) and fll22D for A { *\ eqs. ( g^9|) and ( fOc| ) 



for A-pb and eqs. (|2.45|) for F L . The simple large recoil spin amplitudes in eqs. ([B . 40| ) 
will be used to account for the existence of these relations, allowing us to proceed 
as we did in Sect. 6 of ref. |J. The LO large-recoil expressions are sufficient to 
understand these symmetries, since NLO contributions do not break the pairing of 
the coefficients. 

The following short-hand notation will be used 

2m?, 



M± = 12m^±g 2 (3 + ^), Fee 



.s 



P l = V2Nm B {l-s), P 2 = ^-^ 7=(l-s), 

2v2m^» V s 



and 



d ee Cf M + Sd, C v ee 8C V . (B.45) 



From the simplified expression of the spin amplitudes given earlier in Appendix |B.5.2 
we obtain: 

(2) 

• Ar^' . Being built in terms of just A± and An, the LO behaviour of this observ- 
able can be readily understood (see ref. 0). 

2[C 10 CV + (C 9 + C 7 F)(C 9 , + C r F)]P^e± 



(B.46) 



lr [CI + + (C 9 + C 7 Ff + (C v + C r Ff] P* 4 £\ ' 

where LR stands for "large recoil". Eq. (|B.46|) shows that, at LO, only the 
terms with primed coefficients (SCf) and cross terms like SCiSCf (with i,j = 

(2) 

7,9, 10) might appear in the numerator of Aj, . Neither those involving just 
unprimed coefficients (SCi, with % = 7, 9, 10) nor products of same chirality 
operators (SC^SC^i)) are allowed, whereas in table |7|, the latter are also present 
but they come from NLO corrections and are much smaller than the LO ones 16 . 



Hi 



The fitting coefficients of the terms forbidden at LO are at most 4% of those allowed. 
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Furthermore they are very suppressed by the corresponding NP terms in the 
denominator (see table [§D. Since both primed and unprimed coefficients enter 

(2) 

the full expression of at NLO in the same way, the relations on the first row 
of eqs. ( |2.20|) , ( 2.21 ) and (|2.22|) hold. The remaining relations can be checked 
trivially using eq. (|B.46|). 



dT/dq 2 . The differential decay distribution appears in the denominator of both 
A FB and F L as a sum of 1^ 5Ci 5Cj. At LO it can be expressed as 



dT_ 

dq 2 



LR 



pi 

2q 2 



M, 



(C 9 - C v + (C 7 - C r )FsY + (Cio - CV) P* ? 



+ 



(C 9 - C w + (C 7 - C V )F) 2 + 2(C 9 + C 7 F) (C 9 , + C V F) 
M.{C 2 w + C 2 w )2e A 



2£ 



Although quite cumbersome, eq. (0) allows us to understand table [L(| since 
all the coefficients there appear already at LO. In particular, we can check 
that the largest fitting coefficients (i(o,7) and I( 7)7 )) are enhanced either by the 
square of the factor F = [2m^mB)/q 2 (which becomes very important in the 
low-g 2 region) or by Cg M F, whereas others are enhanced by F but suppressed 
by Cj M (like I(q,9) and /(o,9')) an d the remaining ones are not enhanced at all. 
Eq. (|10D can be used also to verify the relations in eq. (|2.40j) . 



Afb- At LO and in the large recoil limit, the numerator of this observable has 
a structure given by 



1FB 



LR 



-G^P 2 [Ci (C 9 + C 7 F) - C w (Cy + C r F)]g 
dT/dq 2 



(B.47) 



All fitting coefficients in table arise already at LO except for those that 
involve a primed coefficient, i.e. (0, 7'), (0, 10') and (7, 10'). In the case of I(o,r), 
the effect of the enhancement factor F at low-g 2 explained above is particularly 
visible, while I(7,w), which also receives this enhancement, is suppressed by 
Cf M and 

I(o,io') is n °t enhanced at all. 

Regarding eq. (|2.39| ), the first and the second equalities are LO relations due to 
the antisymmetric behaviour of primed and unprimed coefficients in eq. ( p.47|) , 
whereas the last one appears only at NLO but respects the same symmetry. 

Fij. The numerator of this F^ at LO and in the large recoil limit simplifies into 



LR 



2P 2 P 2 [(Cio - CW) 2 + (C 9 - Cy + (CV - C r )Fs) 2 ]g 

dT/dq 2 



(B.48) 
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Using eq. ( p.45|) we can expand the numerator of Fi, into products of NP 
Wilson coefficients. This is enough to derive all relations in eq. ( [2.45|) and to 
explain the enhancement of some fitting coefficients over others in the low-g 2 
region. 
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